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ABSTRACT 

The mean density of the universe is equal to the mass of a large galaxy cluster 
divided by the equivalent co-moving volume in the field from which that mass 
originated. To re-examine the rich cluster $7 value the CNOC Cluster Survey 
has observed 16 high X-ray luminosity clusters in the redshift range 0.17 to 0.55, 
obtaining approximately 2600 velocities in their fields. The systemic redshift, the RMS 
line-of-sight velocity dispersion, ai, and the mean harmonic radius, derived 
for each cluster using algorithms which correct for interlopers in redshift space and 
measure the angular extent of the sampling. The virial mass, and its internal error, 
are derived from these data. The cluster luminosity, corrected to z = 0, is estimated 
from the r band luminosities of the cluster galaxies. Directly adding all the light to 
Mr{0) = —18.5, about 0.2L^,, and uniformly correcting for the light below the limit, 
the average mass-to-light ratio of the clusters is 283 it 27h Mq / Lq and the average 
mass per galaxy is 3.5 it 0.4 x lO^^/i^^ Mq. The clusters are consistent with having a 
universal M^/L value (within the errors of about 20%) independent of their velocity 
dispersion, mean color of their galaxies, blue galaxy content, redshift, or mean interior 
density. Using field galaxies within the same data set, with the same corrections, we 
find that the closure mass-to-light, Pc/j, is 1160 it ISO/iMq/L© and the closure mass 
per galaxy, pc/4){> 0.2L,,), is 13.2 it 1.9 x 10^^/i~^ M©. Under the assumptions that the 
galaxies are distributed like the mass and that the galaxy luminosities and numbers 
are statistically conserved, which these data indirectly support, JIq = 0-20 it 0.04 it 0.09 
where the errors are, respectively, the la internal and an estimate of the la systematic 
error resulting from the luminosity normalization. 
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1. Introduction 



Clusters of galaxies are the largest collapsed objects in the universe and play a particularly 
important role in the problem of estimating the mean density of matter that participates in 
gravitational clustering, po- The standard procedure is as follows. The mass of a cluster, within 
some radius, is measured with any of a number of techniques, such as galaxy kinematics. X-ray 
profiles, or gravitational lensing. To relate cluster masses to po = M/V requires an estimate of 
the co-moving volume, V , from which the clusters collapsed, normally done via the luminosity 
of the cluster galaxies, L, in ratio to the field luminosity density, j, so that V = L/j. The 
value of r^o = Pa/ Pc can therefore be rewritten as the cluster mass-to-light ratio, M/L, divided 
by the closure mass-to-light ratio, {M/L)c = Pc/j- The quantity is independent of Hq, 
and determines the future of the expansion of the universe, in the absence of a cosmological 
constant or a "hot" component of the mass field. The classical cluster mass estimator is the 
virial mass, My, a straightforward global estimator suitable for relatively sparse data (but see 



Bahcall &: Tremaine 198l| ). This subject has a long history, with relatively stable results. That is, 
the virial mass-to-light ratios of clusters, M^/L, are generally in the range 200 — 4OO/1M0/L0 



{e.g. 


Zwicky 19331, Smith 193€, 


Zwicky 1937, Schwarzschild 1954, 


Gunn 1975, 


Ramella, Geller & 


|Huchra 1989|, David, Jones, & Forman 1995, 


Bahcall, Lubin & Dorman 19951) 


which, in ratio to 
icates — 0.2. 


(M/L)c ~ 1500/1 ( 


Efstathiou, Ellis & Peterson 1988, 


Loveday et al. 1993) inc 



The low value of JIq from cluster M^/L measurements is not generally accepted as a definitive 
measurement of the field value of because the random errors are generally fairly large, and 
there are uncontrolled systematic errors in both the mass and the light measurements which 
potentially allow the cluster M^/L ratios to be consistent with Q. = 1. 

There are two latent problems in the M^/L estimate of il. The reliability of the virial 
mass statistic as indicating the total gravitational mass of a cluster is critically dependent on 
whether the galaxy distribution traces the total mass distribution of the cluster. The large color 
differences between cluster and field galaxies means that their recent star formation histories are 
different, which may simply be a recent decline associated with infall into the cluster, or, it may 



indicate drastically different star formation efficiencies at early times ( Dressier 1980 , Dressier &: 
Gunn 19831 , [Butcher fc Oemler 198^ |Yee et al. 1995|) with the consequence that the luminosity 
per unit mass could be quite different between cluster and field. The possibility that the total 
mass distributions of galaxy clusters are more extended than their constituent galaxies has been 



recognized for many years (Limber 1959, West & Richstone 1988, Carlberg 1994). The virial mass 
is a usefully accurate measurement of the mass within the orbits of the galaxies, and completely 
independent of any anisotropy of the velocity ellipsoid, but if the cluster light is more concentrated 
than the cluster mass, the virial mass will be an underestimate of the total mass. This assumption 
of whether the mass and galaxies are similarly distributed is straightforward to test given sufficient 
data to derive radially resolved luminosity and velocity dispersion profiles. 



Dynamical mass measurements have the benefit that their strengths and weaknesses are 
already relatively well understood and are easily studied further with n-body simulations. 
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Although not of direct interest here, dynamical mass estimation can be extended to well beyond 
the virialized region, by measuring peculiar velocities in the infall region, although then the results 
depend on the biasing of galaxies with respect to the mass field. The overall goal of this project 
is to measure the total mass and luminosity contained within the central virialized region of the 
cluster to establish a value of Q, and measure any biases among cluster galaxies, cluster mass 
and field galaxies. Beyond the virialized region there is an infalling mixture of galaxies, gas, and 
dark matter which is very likely statistically identical to the field. Other than the few very old 
central galaxies and a central cD the bulk of the cluster galaxy population is likely to be composed 
of infallen field galaxies. Because our sample of galaxies extends from the cluster core to the 
distant field we can check whether galaxy population modifications only occur in the virialized 
region of the cluster. Therefore we will be able to empirically extend our measurement of from 
the virialized mass distribution to the 17 of the field. The scope of this project lead to a broad 
collaboration under the title of Canadian Network for Observational Cosmology (CNOC). 

The twin complications for interpreting velocities are that cluster galaxies and surrounding 
highly-correlated field galaxies are intermingled in redshift space, and, individually clusters are 
quite nonspherical, and internally clumpy, with the velocities at large radii being influenced 
by surrounding large scale structure. Hence, there are two main requirements for a successful 
measurement: accurate control of the background and a sufficiently large sample to average over 
the aspherical complications of single clusters. A feasibility study ( Carlberg et al. 1994| ) showed 



that approximately 1000 galaxy velocities and a comparable number of field velocities distributed 
over 10 or so cluster fields is sufficient to indicate whether the virial mass is the total mass of the 
clusters. The built-in field sample allows a study of the the relative galaxy populations of cluster 
and field, and the closure mass-to-light and mass per galaxy. 

There is no question that the cluster M^/L technique gives, in principle, an accurate $7 value. 
The key issues are to constrain the technique's systematic errors and minimize its random errors. 
The main areas of concern are the accuracy of the virial mass and galaxy population differences 
between the cluster and the field. Our data allow these problems to to addressed using our 
sample alone, since it contains more than 1000 field galaxies, and offers the considerable benefit of 
homogeneous selection as discussed in Section 2. In Section 3 the cluster RMS velocity dispersions 
and mean harmonic radii are derived. These quantities are used to find virial masses in Section 4. 
The luminosity and numbers of cluster galaxies are used in Section 5 to give M^/L and M^/Nl 
ratios. The field values of these quantities are derived in Section 6, and the resulting Q values in 
Section 7. A companion paper ( ICarlberg, Yee k Ellingson 1995| ) will examine the relative light 
and mass profiles of the clusters. 



2. Sample and Observations 



Clusters at moderate redshifts, z ~ Y3, have a number of advantages for mass estimation. 
They are sufficiently distant that they have a significant redshift interval over which the density 
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of foreground and background galaxies are nearly uniformly sampled in redshift. An accurate 
background estimate is crucial for measurements in the outskirts of the cluster profile. Clusters 
at higher redshifts have increasing fractions of relatively blue galaxies ( [Butcher &: Oemler 198^ ) 
which increases the similarity between cluster galaxies and field galaxies (but they are by no 
means identical), and, because the mean galaxy color varies, the effect of varying galaxy types on 
the total luminosity can be measured. 

There are a number of practical considerations which motivate our choice of clusters. 

At z ~ Vs) 

a comoving Abell diameter of 3h ^ Mpc spans an angle of about 13.3', which 
is sufficiently small that uniformity of photometry and sample selection is relatively easily 
assured. This size is also comparable to the field size of the Canada-France-Hawaii Telescope 
(CFHT) Multiple-Object-Spectrograph (MOS), approximately 10' square ( ILeFevre et al. 1994| ). 



These clusters are also ideal targets for other types of mass measurement observations, such as 
gravitational lensing and X-ray plasma profiles. 

Our cluster sample was selected from the Einstein Medium Sensitivity Survey Catalogue 
(IGioia et al. 199q , [Henry et al. 1992| , pioia k Luppino 1994| ), which is a collection of serendipitous 



objects discovered in X-ray images taken for other purposes. For our survey we chose clusters 
with z > 0.18 (a consequence of the poor blue response of the CCDs then available) and in the 
declination range —15 < 5 < 65, most suitable for CFHT. To select uniformly clusters that are 
likely to be rich in galaxies, have a large velocity dispersion (both of which make the measurement 
relatively easier) and are guaranteed to have a substantial virialized component, we chose those 
with > 4 X 10^^ erg s'^ and /x > 4 x lO'^^erg cm'^ s"^. The A2390 cluster, a rich ([Abelj 
1958), high X-ray luminosity ( Ulmer et al. 19*8^ ) cluster was added to fill an RA gap. 



Observations were made at CFHT in 24 assigned nights in 1993 January, June, October and 



1994 January and March. The observational techniques and data reduction are described in Yee, 



Ellingson h Car lb erg 199^ (hereafter YEC), and the data will be described in a series of papers 



{e.g. Yee et al. 1996| ). The fields (either EW or NS strips between 9' and 45' wide and about 



8' high across the cluster center) were imaged with the Multi-Object Spectrograph ( p^eFevre 



et al. 1994) from which Gunn g and r magnitudes are derived, and masks of spectrograph 
entrance slits were designed for the nonstellar images. The wavelength range of the spectra was 
shortened with band limiting filters chosen to match the cluster's redshift, which typically allowed 
about 100 targets to be observed on a single mask. The limiting magnitude for spectroscopy for 
each cluster is set to optimize the number of cluster (as opposed to field) redshifts obtained. The 
resulting spectra are cross-correlated with a set of templates to give velocities accurate to about 
100 kms~^ in the rest frame of the cluster. The resulting "pie diagrams" for the entire sample 
are shown in Figures || and |2[ where redshift increases with distance from the vertex, and the 
angular co-ordiante is the physical separation in the RA or Dec direction, with the zero co-ordinate 
being at the location of the brightest cluster galaxy (which is not necessarily at the center of the 
observed field). 
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The fraction of galaxies for which we assign shts and obtain redshifts decreases for fainter 
galaxies. Such a selection function, which decreases with magnitude, optimizes the rate of return 
of cluster redshifts, but needs to be carefully measured to reconstruct statistically the properties of 
the complete sample. The algorithm for deciding which galaxies will be assigned spectroscopic slits 
has the highest priority at the expected m* of the cluster. The resulting n{m) of the subsample 
with redshifts is nearly constant with magnitude, whereas the total numbers rise with magnitude. 
The redshift sample is corrected to the photometric catalogue with a magnitude selection weight, 
Wmimr). The limiting magnitude is defined here as the where Wm = 5. The MOS has 
approximately a fixed number of slits per unit sky area (about 100 in 100 square arcminutes), 
hence there is a geometric selection weight, Wg{mr, x,y). For quantities that depend on surface 
area, these weights need to be supplemented with the fraction of circle sampled by galaxies at 
varying radii from the chosen centre. The defining relations for these weights are described in 
YEC. All the results given below incorporate these weights as is appropriate. 

In this paper, all distance dependent quantities are calculated assuming 
Hq = lOOkms"^ Mpc~^ and qq = 0.1, A = 0. The choice of qo is for approximate consistency with 
the most straightforward interpretation of our results. The dependence of distance and luminosity 
on go largely cancels for the evaluation of O because it is defined relative to the field in the same 
redshift range. Our luminosities are corrected to a redshift zero for the sake of approximate 
comparability with low redshift results. These corrected quantities should only be compared 
to others that our adjusted to our photometric system. Hence, care should be exercised in any 
direct comparison of our corrected z quantitites, for instance mass-to-light ratios, to similar 

quantities at z = 0. 



3. Cluster Dynamical Parameters 

The redshifts and positions of galaxies are used to define a characteristic velocity and a 
characteristic length scale of the cluster. As a consequence of the virial theorem, the RMS 
velocity dispersion has the important property that in a spherical system its value is completely 
independent of any variation in the shape of the velocity ellipsoid. For a triaxial object the 
line-of-sight RMS velocity does depend on viewing angle, which is precisely why we want to 
average over a dozen or so clusters. The line of sight velocity dispersion of a cluster is defined as 

^?=(e^') (1) 

where the Avi = c{zi — 'z)/{l + 'z) are the peculiar velocities in the frame of the cluster and z is 
the weighted mean redshift of the cluster. The weights used in Eq. ^ are magnitude dependent 
geometric weights. Unweighted velocities give similar results. 
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The virial mass estimator normally uses the projected mean harmonic pointwise separation 



where the ij sum is over all pairs. Being a pairwise quantity, Rh is sensitive to close pairs and is 
quite noisy (Bahcall & Tremaine 1981). Furthermore there is no straightforward way to correct 



Rfj for the approximately rectangular "window" which encloses our cluster sample. That is, 
simply weighting with the fraction of a circular aperture that is enclosed within our rectangle 
leads to a systematic underestimate of Rh- 

Here we introduce an alternate estimate of Rh- A pair of galaxies at projected co-ordinates rj 
and Tj are statistical representatives of all galaxies at those projected radii. Under the assumption 
of axial symmetry the angle between these two vectors is a uniform random variable. Hence, 
some immediate averaging is possible by imagining that one (or the other) of the particles has its 
mass distributed like a ring, with a radius equal to its radial location with respect to some cluster 
center. More formally, the expectation value for random angles of the pairwise potential l/jr; — rj] 
is the potential between a point and a ring, which is (noting that the result is independent of the 
sign of the cosine function for an integral over a circle) 



RK' = (e-^) T^^^^^^fJ- 



de 



\-2 2 



(3) 



where kf- = Arirj/{ri + rj)'^ and K{k) is the complete elliptic integral of the first kind in Legendre's 



notation ( Press et al. 1992 ). The quantity Rh will be referred to as the ringwise projected 
harmonic mean radius. Unlike the original Rh this modified Rh requires an explicit choice of 
the cluster center and assumes that the cluster is symmetric about the center. Neither of these 
pose any practical difficulties. The profile analysis will make these assumptions, but it should be 
noted that for flattened clusters the resulting Rh will be a small overestimate if the field sampled 
happens to lie along the major axis of the cluster. 

There are two substantial benefits to be had from the definition of Rh in Equation ^. For 
close pairs the divergence is logarithmic instead of 1/r, which makes Rh less noisy than Rh- Of 
immediate practical interest here is that the value of Rh is readily determined for data sets where 
a strip (sometimes of varying width) across the center has been sampled. For instance, any of our 
strips can be artificially narrowed to test how much Rh varies. For a factor of two reduction in 
the width of the A2390 data, Rh declines from 348 (already a substantial underestimate) to 267 
arcseconds, whereas Rh drops from 533 to 517 arcseconds. 

The three dimensional virial radius, r^, is ry = 'KRh/2 (iLimber Mathews 196^ ). With cji 
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and the virial mass of a cluster is calculated as 



= ^(^fr^- (4) 



3.1. Cluster Membership 

The complication in applying Equations |l] and ^ to redshift data is to decide which galaxies 
belong to the cluster and which are likely to be interlopers in redshift space. Because galaxy groups 
and other clusters have such a high clustering probability near a rich cluster this background is 
not expected to be smooth. Our magnitude distribution of galaxies with redshifts was designed to 
maximize the abundance of cluster galaxies and minimize the background. That is, the number 
of field galaxies per magnitude, n(m), rises steeply with magnitude, n{m) oc 0.4m, whereas the 
cluster luminosity function is much shallower, n{m) oc 0.25m or less, below m*, the characteristic 
brightness of a cluster galaxy. For each cluster we estimated and then set our spectroscopic 
exposures so that we expected to obtain redshifts down to the magnitude where cluster galaxies 
are starting to become less numerous than field galaxies, which typically occurs 2.5 magnitudes 
below m^,. The higher priorities are assigned to galaxies with brightnesses near m*. To meet 
the goals of this project it is far better to have the same total number of cluster redshifts in 
many clusters, rather than in a single one. Besides being an efficient procedure, this strategy also 
ensures that the cluster has the maximum possible contrast in redshift space from the field, as 
shown in Figures || and ^. A sampling which is more complete or deeper would mainly provide 
more field galaxies. If more velocities are desired it is far better to increase the sample of clusters 
(to increase the statistical averaging over cluster to cluster variations, such as substructure) rather 
than decreasing the random error in a few clusters. However, fewer than 50 velocities in a cluster 
often leads to large uncertainties because the velocity structure of the cluster is not clear, in which 
case one tends to overestimate the velocity dispersion. 

The calculation of the RMS velocity dispersion of a cluster uses only the redshifts with no 
positional information. Defining the extent of rich clusters in redshift space is a relatively difficult 
problem in general {e.g. Bird Sz Beers 1993 , Bird 1995| ) because cluster members moving at a 



large velocity with respect to the center of the cluster are impossible to distinguish from field 
galaxies. Furthermore, many statistical tests to assess the membership status of galaxies in the 
tails of the velocity distribution are based on the assumption that the underlying distribution is 
a unimodal Gaussian, which is generally not true for such a kinematically complex system as an 
accreting cluster and the surrounding large scale structure. We have reduced the problem through 
the magnitude distribution of our cluster sample, which provides data that leaves clusters as 
cleanly defined as possible in redshift space. 

The procedure for defining the redshift range of the cluster is a manually iterated procedure, 
roughly equivalent to measuring the width of an emission line in a spectrum. First, an initial 
estimate of the outer limits of the redshift range of the cluster is made, using Figures ||, |^, and 
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^. The redshift limits adopted are shown as dotted lines in Figure |3[ The initial estimate of the 
redshift range is used to calculate a trial velocity dispersion, ai. A second algorithm is used to 
assess the validity of this velocity dispersion in the presence of interlopers in the outer parts of 
the redshift space of the cluster. The velocity dispersion validation algorithm works as follows. 
The weighted (as in Eq. |T|) velocities within 15c7i (or the edge of the band-limiting filter if that is 
smaller) of the cluster center are put into bins O.lfii wide (within reason the bin width has no 
effect). The mean density of the background in velocity space is calculated from the mean bin 
density > 5(Ti away from the cluster. This background is then subtracted from the data within 3ai 
of the cluster center and the velocity dispersion a[ is calculated, along with a Bootstrap estimate 
of the error in a'l- If a'l and cJi are within one standard deviation, ui, is accepted as the cluster's 
velocity dispersion; if they are not equal, the redshift limits are further adjusted. 

The adopted ranges in redshift, as defined by the above procedure, and radius, as defined 
by the rectangular region which bounds the data, are given in Table |l[ Column 2 gives the 
weighted mean redshift, columns 3 and 4 gives the dimensions of the bounding box in RA and 
Dec as converted to physical lengths, Dj^{z)A9, at the redshift of the cluster. Columns 5 and 6 
give the lower and upper redshift boundaries. In Table ^ the second column gives the ratio of 
the "unlimited" a'l to the value accepted, and the third column gives the 1 standard deviation 
confidence interval for the cr'i/ai ratio. Each cluster's velocity dispersion is calculated using 
Equation || with the weighted data in the cluster's redshift range given in Table |l[ The results are 
given in Table Column 2 of Table ^ is the redshift, column 3 is the number of galaxies above 
the magnitude limit in the cluster RA, Dec and redshift range. Column 6 is ci. The sky positions 
of the cluster galaxies for the inner 1000 arcsec of the sample are shown in Figure ^. 

The obvious redshift limits were acceptable in about half of the clusters. In the other half 
the limits were adjusted until a statistically acceptable result was found, or until no stable 
value emerged. The velocity dispersion checking algorithm indicates (Table ^) that the velocity 
dispersion of MS0906+11 is significantly too large. No consistent solution could be found for 
this cluster, which appears to be an indistinct binary in redshift space. All other clusters have 
velocity dispersions that are well within two standard deviations of the input values. Based on 
this test, we will accept the results of Table ^ as our standard values. One could apply the slight 
corrections of Table ^, but that is simply changing the values within their standard error. The 
cluster MS1358+62 appears to be an unequal binary in declination (Fig. |2|), although its velocity 
dispersion appears to be acceptable. 

The modified mean harmonic radii, Rh, are calculated for each cluster using the entire angular 
extent of the available data in the redshift range of the cluster as fixed for the ai calculation. The 
mean interior overdensity inside r„ is an indicator of the radial range of the cluster that has been 
observed. The angular extent of the clusters is defined by the angular limits on the sky of our 
sample. 



It is expected (Gott fc Gunn 1972) and borne out by numerical simulation (e.g. Crone 
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Evrard & Richstone 1994) that, on the average, substantiahy all the virialized cluster mass is 
contained inside the radius where p > 200pc- Outside that radius the velocities very rapidly 
become dominated by infall. Therefore the ratio p(r^)/pc is used as a test of the completeness of 
the radial extent of the sampling of the cluster. This ratio is calculated as the mean density inside 
Tt, to the critical density at the observed redshift by dividing the virial mass (Eq. ^ by the volume 
inside r^. 



p(r^) _ 1 3My 

(5) 



Pc{z) pc{z) Aixrl 

6 al 



{l + zY{l + QQz)Hlrl- 

The critical density at redshift z is Pc{z) = pc(0)(l + z)^(l + Q.qz), which for = 1, has the 
(1 + zY dependence of the mean density. For < \ Eq. ^ means that the density ratio at z 
is somewhat larger than for fig = 1, with the difference being 25% at z = Ya for our adopted 
2go = 1^0 = 0.2. 

The total virialized mass of the cluster should be accurately estimated if the mean density 
inside r^, is 200/3c or less. Table ^ includes 6 clusters for which this criterion is met. For reference, 
the Coma cluster has Rh ^ 1.78/i~^ Mpc (1.48°, [Schwarzschild 1954| ) and, with cti of 1040 kms~^ 



(Peebles 1970D has p(r^) = 77pc- 



3.2. Error Analysis 

The errors in ui, Rh^ M^, and M^/L are assessed using the Jacknife technique. This is one 
of the simplest resampling techniques, wherein partial standard deviations, 6i, are calculated by 
taking the difference between the /(xi, . . . where / is any statistical quantity calculated 
from the entire data set, and the same quantity calculated dropping one element of the data set, 
5i = f{xi, . . . , Xn) — f{xi, . . . , Xj-i, . . . , Xn)- The estimate of the variance is Jn/ {n — 1) J2i 



( Efron 1981 , Efron fc Tibshirani 1986 ). For a Gaussian distributed sample the Jacknife variance 
converges to the normal value. A closely related technique is the Bootstrap, in which randomly 
drawing with replacement from the original data set creates many new data sets of the same size 
as the original which are then analyzed in precisely the same manner as the original to give a 
distribution of results from which confidence intervals are calculated. A significant difference is 
that the Bootstrap data sets have repeated values, which for pairwise statistics gives singularities. 
Also, the Bootstrap can be computationally quite expensive, since tests indicate ( [Efron fc 



Tibshirani 1986) that 300 or more resamplings are required to give solid results. The errors given 
in Table |3| (in the columns labeled and e^-) are approximately in accord with a Vn expectation, 
but have significant differences in detail. The reliability of these errors is extremely important 
for one of the major results of this paper. We believe the pragmatic nature of the Jacknife error 
estimate is well suited to the task. 
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3.3. Sample Commentary 

A cluster-by-cluster comparison with the results of other workers is done in the data papers 
for each cluster. In a number of cases our velocity dispersions are notably lower than those 
previously found. For instance, we find for Abell 2390 that ai = 1100 it 63kms~^ whereas a 
previous study in the central region found 21121^q7 kms^^( [Le Borgne et al. 1991 ). In general 



the differences from previous work can be attributed to three factors. First, the precision of 
our velocities, ~100 kms~^, reduces catastrophic velocity errors and helps identify significant 
substructures. Second, our data cover a large radial range in the cluster, which makes them less 
subject to local substructure, and more representative of the RMS value. Third, the bluer galaxies, 
which often contain measurable emission lines, statistically are found to have a higher velocity 
dispersion than the redder absorption line galaxies, an effect which is particularly prominent near 
the projected center of the cluster (Abraham et al. 1995, [Carlberg, Yee &: Ellingson 1995| ). Since 



redshifts are easily measured from emission lines this can lead to an upward bias in the velocity 



dispersion, but which is demonstrated to be quite small in our data ( |Yee, Ellingson fc Carlberg 



19961) . These effects can be deadly in the presence of substructure, for which there is evidence 
in many of our clusters. In general the clusters are quite regular, possibly connected to their 
high X-ray luminosity, and usually does not upset any of the results below. The exceptions are 
identified and removed, as discussed below. 

The quantity p(r„) in Table |3| varies from a low of 84pc in MS1621+26 to 1135pc in 
MS1455+22, discounting MS0906+11. It is clear that 'p{rv) is high for those clusters for which 
only a central field was observed. In particular, we expect on the basis of Coma data and n-body 
simulations, that it is necessary (but not sufficient) that p(r^) should be less than 200pc if the 
formal (Eq. ^ is interpreted as the total virialized mass. li'p{ry) > 200pc then it indicates that 
r„ is an underestimate of the true value for the entire virialized system, which is a straightforward 
consequence of the limited angular sampling of some of our clusters. This undersampling is an 
important issue, but is not a significant problem for the major goals of this paper because some 
quantities, such as My/L, are insensitive to the outer limits of the sampling, and it is possible 
to define a characteristic radius at a fixed mean interior density, which is insensitive, within 
reasonable limits, to these sampling variations. 



4. Virial Mass Analysis 

The quantities ai and r„ are given in Table ^. The resulting My (Eq. ^) are given in Table |^, 
where column 2 repeats the 'p{ry) column of Table |3|, column 3 is My in units of solar masses, 
column 4 is the total observed luminosity in units of solar luminosities (see Section 5), columns 5 
and 6 are the canonical My/L ratio and its standard error from the Jacknife technique. It must 
be borne in mind that the masses given are defined by the radial extent of the sample, and must 
be referred to the mean density given in column 2. That is, if the cluster extends beyond the size 
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of our observed field, our mass is an underestimate of the total mass. Nevertheless, it remains of 
considerable interest to compute an M^jL for the part of the cluster we observe. Furthermore it 
is relatively straightforward to bring all the masses to a common ground. 

The virial mass, as calculated from a data set which misses the outer parts of the cluster, 
will still accurately estimate the mass contained within the orbits of the galaxies observed, with a 
weak dependence on the radial extent of the sampling as we show here. The virial theorem is a 
derived by multiplying the Jeans equation of stellar "hydrostatic" equilibrium (a vector equation) 
with co-ordinate vectors and integrating over the volume of the system {^e.g. [Binney &: Tremaine) 



19871) . The outer limit of the volume integral is usually taken at infinity, but, when evaluated for 
a finite subset of the entire volume, gives an additional term to be evaluated at the outer surface. 
The scalar virial theorem for a spherically symmetric system bounded at radius becomes 

''\al + 2al)pdV + r Z^f^^pdV = inp^rl (6) 
Jo r 

The right hand side of Eq. |6| is sometimes called the 3PV surface term. The effect of ignoring the 
surface term for a gravitating system is to overestimate the mass of the system, since the surface 
pressure reduces the amount of mass to keep the system in equilibrium. 

A singular isothermal sphere gives an estimate of the largest reasonable error if we use Eq. ^ 
on a partially sampled cluster. The isothermal sphere provides an upper limit, since equilibrium 
stellar systems normally have a falling velocity dispersion with radius, so the isothermal sphere 
has an unusually large surface term. The density profile of a singular, isotropic, isothermal sphere 
of velocity dispersion a = cr^ is p = (27rG)~^cr^r~^. Equation |6|then becomes 

3^2^^ _ GMl ^ ^2^^^ 
n 

where Mb is the mass inside the boundary. If the surface term is neglected, then the derived 
virial mass will be an overestimate of the mass contained inside r^, SG^^a'^rh, instead of the 
correct 2G~^a'^rh. If the density profile declines more rapidly than r~^, as is indeed the case 
for the clusters studied here, then the error introduced through the neglect of the surface term 
is reduced. The isothermal sphere has the pleasing property that the relative size of the surface 
term is completely independent of shape or size of the surface. We conclude that if only an inner 
part of the cluster profile is sampled, then the calculation of the virial mass will be at most a 50% 
overestimate. 

Under the assumption that clusters are approximately p{r) oc r^^ (or any other mass profile 
of interest), the masses can then be cautiously extrapolated to other radii or mean interior 
densities. The results of extrapolating the masses to a constant mean interior density of 200pc 
using M200 = My y^'p(ry)/{200pc) are given in the second last column of Table 

In "raw" M^ the most massive cluster in the sample is A2390 (discounting MS0906-I-11), with 
My = 1.7 X lO^^h^^ Mq, which makes it comparable to Coma's 2.1 x lO^^h"^ Mq. Coma and 
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A2390 have similar mean interior overdensities at r„, 77 and 110, respectively, putting the mass 
comparison on nearly the same density basis. The masses of the other clusters range a factor of 
about 2.5 bigger and smaller than M200 ~ 5 x 10^^ Mq. This rather narrow range in mass is a 
consequence of our X-ray luminosity and flux selection. It is perhaps not too surprising that one 
of the most massive clusters, MS0016+16, is also the highest redshift one in our sample. 



5. Cluster Galaxy Luminosities and Numbers 

There are two practical approaches to relate empirically the cluster contents to a co-moving 
volume in the field. A third approach, a total cluster baryon inventory ( [White et al. 1993| ), cannot 
at present directly relate the cluster to the field since most field baryons are in some unobserved 



form and one must turn to a Big Bang Nucleosynthesis argument ( Walker et al. 1991 ) to obtain a 
field density of the baryons. However, the large reservoir of non-stellar baryons in clusters raises 
concern about the assumption that the fractional conversion of gas into stars has been the same in 
clusters and the field. Our galaxy sample extends from cluster to the far field, allowing a number 
of tests for variation. The two "binary" clusters, MS0906+11 and MS1358+62 are excluded from 
all averages calculated below. 

The total luminosity of a cluster estimates the co-moving volume from which it collapsed, 

V = L/ j. This approach assumes that the total luminosity is conserved, or that any systematic 
change in galaxy luminosities during the gravitational assembly of the cluster can be measured 
and corrected. However the subject of galaxy evolution within clusters remains quite controversial, 
at the very least because the observational phenomena are not yet well defined. 

An alternate approach to measuring the equivalent field co-moving volume of the cluster 
is to use the numbers of cluster galaxies above some absolute magnitude in ratio to the field, 

V = pc/4'{> L). Each of these methods has different strengths and weaknesses, but they are 
substantially independent since cluster total luminosities are dominated by the brightest galaxies 
but the total numbers are dominated by the low luminosity galaxies. One possible evolutionary 
effect is the merging of the stellar components of two galaxies which makes no difference to 
the total luminosity if there is no accompanying episode of star formation. Merging does, of 
course, decrease the numbers of galaxies. It is well known that cluster galaxies are on the average 
significantly redder than field galaxies (Figure ^), which is taken to indicate that star formation 
has decreased in cluster galaxies and hence they are less luminous per unit stellar mass than 

in the field. Fading will decrease the numbers of cluster galaxies above some fixed luminosity 
limit relative to the field, although, if the limit is well below M^,, then the decreased numbers 
will be a small fraction of the change in luminosity. For instance, if all cluster galaxies down to 
Mj.(0) = —18.5 were included in the total count, but there had been a one magnitude fading with 
respect to the field, then galaxies down to —17.5 should have been included in the total count. 
The expected increase in numbers (for a Schechter function with = —20.3 and a = 1.25) would 
be only 17%, as compared to a 150% correction required to the luminosity. Cluster galaxies are 
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not simply post-star formation field galaxies with diminished luminosities, since clusters contain a 
population of very luminous, but very red. Elliptical and cD galaxies in their centers. These may 
have been created through mergers, or they could have tapped the large baryon reservoir present in 
the cluster to create more stars. If galaxy luminosities are either increased or decreased, the L of a 
cluster responds in proportion, but the numbers of cluster galaxies above some magnitude change 
far less. The evolution of cluster galaxies relative to field galaxies, especially at these redshifts, is 
sufficiently ill defined that we do not, at this stage, want to make a highly uncertain differential 
correction. However, we will use mass-to-light and mass-per-galaxy to obtain two substantially 
independent estimates of Q, which in turn will be used to address the notoriously difficult problem 
of constraining the size of systematic errors. That is, if the two estimators give discrepant values 
beyond the errors, then we would have detected a systematic effect (whose source would then be 
identified and corrected, if possible). We cannot empirically detect systematic errors smaller than 
our random errors which then becomes our systematic error estimate. 

The observed luminosity of the cluster galaxies, summed to Mr{0) = —18.5, corrected as 
described below, is reported in Table ^. The luminosity is measured in the Gunn r band, using 
Mr{Q) = 4.83. We also note that r band selection has the useful benefit that it is not unduly 
sensitive to star formation increases, although these are readily detected in the g — r colors. The 
measured are K corrected using the following method. Model K-corrections in r and g — r 
colors as a function of redshift for non-evolving galaxies of 4 spectral types (E-I-SO, Sbc, Scd, 
and Im) are derived by convolving filter response functions with spectral energy distributions 
from [Coleman, Wu &: Weedman (1980) . These values are then corrected from the AB system 



to the standard Gunn system ( Thuan &: Gunn 1976| ). For each galaxy with redshift, a spectral 



classification is estimated by comparing the observed g — r color with the model colors at the same 
redshift. The spectral classification, obtained via intepolation, is treated as a continuous variable 
within the 4 spectral types. From the spectral classification, the appropriate K-correction to the 
r mangitude is then derived using the models. Our data strongly indicates that the K corrected 



cluster luminosity function is brightening with redshift ( Yee, Ellingson fc Carlberg 1996 ), hence 
our luminosities are given a mild evolutionary correction, E{z) = —z ( Yee &: Green 1984 ), so that 
they are appropriately faded to the current epoch. 

The corrected magnitudes and colors are shown in Figure ^. An entirely empirical redshift 
normalization of the colors, {g — r)z = {g — r)/(1.2 -|- 2.33 * {z — 0.3)), is shown. This relation is 
based on a linear fit to the redshift dependence of the red sequence of the bright cluster galaxies. 
Because the median color of these clusters is always quite red, there is little range in the quantity 
{g — r)z- The fraction of the galaxies that are blue can be estimated with a parameter similar, 
but quite different in detail, to the Butcher-Oemler blue fraction ( [Butcher fc Oemler 198^ ). The 
quantity Fb here is the fractional weight of all galaxies bluer than (g — r)^ = 0.7. Note that this 
includes all galaxies used in the total luminosity and does not restrict the measurement of the 
blue fraction to some fiducial region. A more detailed analysis of the Butcher-Oemler effect will 



be presented elsewhere (a preliminary account is in Yee et al. 1995 ) 
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It is conventional to quote M^/L values which are based on the "total light", that is, a 
correction is made for the part of the luminosity function below the limiting magnitude. This is 
not strictly necessary here, since the field sample is treated in precisely the same manner, but 
for the sake of comparability we do this correction, for both the cluster and field luminosities. 
The "unobserved light" beyond the limiting magnitude our sample is estimated assuming 
M=K = —20.3 + 51og/i and a = —1.25 in a Schechter luminosity function. A useful limiting 
magnitude is about = —18.5, which is about 0.2L=k. The figures and tables of this paper use 
this limiting luminosity, unless otherwise stated, which helps put all the M^/L values on the same 
basis. There is very little difference in the estimated L's if we include all the galaxies to the 
sampling limit. The estimate from the M^/L has no dependence on the correction to the total 
luminosity, since the field galaxies will have precisely the same multiplicative correction included. 

The luminosities, like the masses, are the total luminosities for the sample, and must be 
referred to its mean interior density. The luminosities can be extrapolated to a common mean 

density, -L200! i^i the same manner as the masses. That is, L200 = LyJ'p{r^) / {2{}Qpc) which builds 

1/3 

in a correlation between -L200 and , but not with the dynamical range observed here. The 
cluster quantities are defined on the basis of a constant ratio with respect to the critical density 
at the redshift of the cluster, Pc{z)^ which increases with redshift. Singe cj^ oc /?^/^M^/^, with the 
redshift dependence of p included this becomes cii oc (1 + z)^/^(l + Q.qzY/^M'^/^ ^ Figure ^ plots 
the redshift adjusted luminosity against the velocity dispersion. This plot excludes MS0906+11 
and MS1358+62. 

5.1. Cluster Mass-to-Light Ratios 

It is immediately apparent from Table ^ that the M^/L do not show much variation outside 
their errors, in spite of the variations in the sampling radii. The mean M^/L for the sample 
(excluding MS0906+11 and MS1358+62) is 283 ± 27/iM0/L0 (in Gunn r, corrected to z = as 
described above), and 245 ± 32/iM0/L0 for the six clusters that have a mean interior density 
less than 200pc) excluding the asymmetric MS1358+62 cluster. The M^/L values as a function of 
redshift are shown in Figure ^. It should be noted that luminosity evolution has been included, 
at the rate of AM = —z. In the absence of this correction the luminosities at z = 0.55 would be 
1.66 times higher, and would lead to a noticeable gradient in the diagram. The corrections to the 
luminosities are all multiplicative factors that are applied both to the field and the cluster data, 
so they have no effect on the Vt value derived. 

The distribution of M^/L differences from the 14 cluster sample mean, normalized to their 
measurement errors, is also displayed in Figure |^ The cluster at —3.6 standard deviations, 
MS1621+26, is one of the clusters with the best data and is also one of the bluest clusters, 
possibly indicates that its luminosity is higher than it would be if the galaxies were allowed 
to age enough that the colors became comparable to other clusters in this redshift range. The 
redshift- Mt,/L plot of Figure ^ also gives the impression that all the clusters have the same M^/L 
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(removing MS0906+11 and MS1358+62) within the errors. The per degree-of-freedom is 1.8, 
but removing MS1621+26 as well reduces it to 1.03 which is consistent with no variation beyond 
the errors. That is, the M^/L of Table ^ are consistent with a universal underlying M^/L, after 
a 3 standard deviation clipping is applied. This means that the measured cluster light is a very 
good indicator of the mass contained within the orbits of the galaxies. It is intriguing to note that 
this appears to be true over a substantial range of mean interior densities which characterize the 



cluster samples. A similar result has been found for X-ray measurement of cluster masses ( pavid 



Jones, fc: Forman 1995| ), although their M/L profiles values decline with increasing overdensity, 
whereas our integrated values are nearly constant, or rise slightly, with mean interior overdensity, 
within the limits expected from the virial surface term of Eq. |^. 

The percentage errors in the M^/L values are in the 10-30% range, depending on the number 
of cluster galaxies. With more data it seems quite likely that it will become evident that there are 
intrinsic variations in cluster My/L values. In particular the cluster light could have a "second 
order" luminosity correction which would depend on the precise mix of galaxies in the cluster. 
This is already the case with MS1621-I-26 which has 98 members with redshifts and is somewhat 
"blue" compared to the trend of colors with redshift. Of course any true M^/L variations cannot 
be bigger than the 20-30% error bounds found here. 



5.2. Cluster Mass-to-Number Ratios 

The mass per galaxy for all galaxies brighter than some chosen absolute magnitude is less 
sensitive to differential fading (or possible brightening) between the field and cluster than the mass 
per unit luminosity. Merging could alter the galaxy numbers, but, in as much as mergers lead to 
ellipticals (the E-l-A galaxies appear to be mostly disk systems. Dressier et al. 1994]) , it does not 



drastically alter the numbers of galaxies. Both fading and merging will cause the M^/Nl ratio 
to overestimate the value of fi. The individual mass-to- number ratios to a limiting magnitude of 
Mr{0) = —18.5 are shown in Figure ^ and given to —18.5 and —19.5 in Table ^. The ^or a limit 
of —18.5 is very large, about 3 per degree of freedom. At a limit of —19.5 the per degree of 
freedom drops to 0.6 if MS1621-I- 26 is dropped. 

Like My/L, My/N shows no significant dependence on other cluster parameters. Figure ^ 
weakly suggests an increase in scatter as a function of the mean interior overdensity, although this 
is likely to be primarily a statistical effect, since the highest density clusters are those that are 
least well sampled. There is a tendency for the high velocity dispersion clusters to have the largest 
M/Nl values. Although this trend is not very significant, it is opposite to a standard biasing 
scheme, "peak-background", in which galaxy formation is initiated earlier and is more advanced 
in clusters of increasing mass. 

The average values of M/N^ for all the clusters are given for a range of limiting magnitudes 
in Table ^. The table demonstrates that the M/N^ is not very sensitive to the limiting absolute 



- 16 - 



magnitude if it is at least 2 magnitudes less than M*, but that the sensitivity rises quickly 
(exponentially is expected for a Schechter luminosity function) near M^. The reduced of 
the deviations from the mean are very large for low luminosity limits. In the case of a limit 
of Mr{0) = —19.5 the reduced drops to 0.8. This is partially a consequence of the smaller 
number of cluster galaxies, typically 60% of the full sample, leading to an increase in the random 
errors. None the less, the reduction in is so striking that it suggests that the most luminous 
cluster galaxies have far less cluster-to-cluster variation than the lower luminosity ones, which is a 
comment that has been made in other contexts. 



6. Light and Number Densities of Field Galaxies 

The closure mass-to-light ratio is defined as {M/L)c = Pc/j, where pc = Sifg /(SvrG) , and 
j = L(p{L) dL, where (p{L) is the field luminosity function. The integration is carried to 
low luminosities by assuming a Schechter form for the luminosity function for galaxies below 
the absolute magnitude cutoff. Similarly the closure {M/Nl)c = pc/4>{'> L), where Nl is the 
number of galaxies more luminous than L. At low redshift {M/L)c is estimated to be about 
15OO;|[j{]/iM0/L0 ( [Efstathiou, Ellis fc Peterson 1988| , [Loveday et al. 1993| ). This value is based 



on a blue selected sample, and integrates all the light in the luminosity function to very faint 
magnitudes. Our selection is based on Gunn r, and we do not sample galaxies much fainter than 
Mr — —18.5. These differences would lead to substantial systematic errors in $7 if we tried to ratio 
moderate redshift cluster light to a low redshift field estimate of the closure density. 

Our survey includes a built-in field sample which can be used to estimate accurately (M/L)c 
and (M/N)c, on the same manner as our cluster values. Each field galaxy is K corrected and 
evolution corrected, and within a redshift range we correct the total luminosity for light below 
our magnitude limit in precisely the same as the cluster galaxies. Because all of these factors are 
multiplicative, they cancel in the eventual ratio of interest for fi. The corrections are done for 
the sake of providing numerical values having some degree of comparability to true low redshift 
quantities. The luminosity density of the field can be nearly trivially calculated from our field 
sample, and offers the benefits that it is selected in precisely the same way, in the same filter band, 
in the same redshift range, and to the same limiting magnitude as the cluster data. This approach, 
using a self-consistent luminosity system for the cluster and field, greatly reduces concerns about 
systematic differences between the low redshift estimate and our medium redshift data. 

Table gives the results of the {M/L)c and {M/Nl)c calculation to the indicated limiting 
magnitudes, averaged over the redshift range 0.2 < z < 0.6. The cluster redshift ranges of Table I 
are widened by Az = 0.01 at both the upper and lower redshift to eliminate any concerns about 
an overestimate of j from a higher density of galaxies in the vicinity of these very rich clusters. In 
fact there is no significant change in the field luminosity density without this extra cut around the 
clusters. The weighted luminosities of the field galaxies, corrected to z = precisely as the cluster 
data, are simply summed in each redshift bin, cutting the sum off at various absolute magnitudes. 
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The volumes are in co-moving h Mpc integrated over the accessible redshift range of the bin 
for our chosen qq, times the total solid angle of our survey, 2367 square arcminutes. There is no 
statistically significant redshift trend in either the luminosity or number density in our evolution 
corrected magnitude system. If evolution correction were not included the values would be about 
35% smaller, and an evolutionary trend would be marginally significant in the luminosity density. 

The random errors of our closure {M/L)c (column 4 of Table |^) are calculated using the 
Jacknife technique assuming that each cluster field should be treated as a data element. This might 
be a conservative estimate of the errors, but fully accounts for large scale structure variations. 

The colors of the field galaxies are shown in Figure |5|. The field galaxies are always 
substantially bluer than the cluster galaxies: the average {g — r)z is 0.55, with no significant 
change with redshift over the entire range of 0.1 < z < 0.6. The cluster galaxies have average 
(5 ~ f)z values ranging from 0.72 to 0.94. Although there is some overlap of the populations, there 
is no overlap of the averages. Hence it is not possible to empirically estimate from this sample the 
luminosities of the cluster galaxies if they had the same colors as the field galaxies. 

7. n 

The value of is estimated under the assumption that the M^/L of the clusters is the same as 
M/L in the field. Our best estimate, in the sense of smallest random errors is the M^/L estimate 
to a limiting magnitude of —18.5 for which we find that = 0.25 it 0.05 at a mean redshift of 
0.31. Brighter absolute magnitude limits give consistent values of fi, but with increased errors. 

The value of 0, from the My/N^ calculation is 0.28 it 0.08, for absolute magnitude limits 
of both Mr{0) = —18.5 and —19.5. This value is consistent with the one found from the M^/L 
argument. We will prefer the M^/L value of mainly because of its smaller errors, but also 
as the conventional approach. The fact that the M^/L and the M^/Nl routes to calculating Vt 
gives similar results is considerable grounds for confidence in the result. The two O values have 
substantially different dependencies on the evolution of cluster galaxies, and are weighted to 
different parts of the luminosity function. Certainly we would not be able to detect any systematic 
errors in cluster galaxy evolution which is smaller than the quadrature summed random errors 
of the two Vl results. We will use this error sum as our estimate of the likely size of any residual 
systematic errors in our cluster Q values, due to galaxy evolution alone. It does not measure any 
systematic dynamical errors in the masses. 

Tables ^ and ^ both have a a common trend M/L values where fainter imposed cutoffs, which 
use more of the real data, find less light than the adopted luminosity function does, although the 
trend is not very significant. 

Since JIq = ^z/{^ + z{l — fi^)), for our average redshift of 0.31 this is VLq = 0.20 it 0.04, 
where the formal 1 standard deviation random error is given. The strength of this value is that 
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is entirely from our survey alone, with its well understood sampling statistics, requires no low 
redshift comparisons, and is done entirely within one photometric system which requires only 
small redshift corrections. The main issue which is addressed in the companion paper ( Parlberg] 
Yee &: Ellingson 1995| ) is whether the virial mass correctly estimates the total mass of clusters, 
with further examination of the differences between cluster and field light per unit mass. 



8. Conclusions 

The CNOC cluster sample is a uniform. X-ray selected, sample of clusters having 
0.17 < z < 0.55. We derive velocity dispersions and a characteristic radius for each of these 
clusters, and calculate their virial mass, based upon objective tests to estimate the radial and 
redshift extent of the cluster members. The Gunn r band luminosity of the clusters is corrected to 
z=0, and the ratio M^/L to a limiting magnitude of —18.5 (corrected to a total luminosity with a 
Schechter luminosity function), is found to have a mean value of 283 it 27/iM0/Lq, comparable 
to low redshift clusters. The clusters are consistent with having identical My/L values (for the 
corrected luminosities) within the sample variance of ztlOOMQ/L©. This is consistent with 
clusters being approximately isothermal, and being dominated by a relatively red population of 
galaxies at all the redshifts observed. 

The field luminosity density is calculated from the same data set so that we obtain an {M/L)c 
in the same redshift range. We find {M/L)c = 1160 it ISO/^Mq/Lq, where the error is estimated 
from field-to-field variations. This is consistent, within the errors, with low redshift results, but it 
should not be compared directly to them because the two are not derived in an identical manner. 
We conclude that the cluster virial masses universally indicate that Vt = 0.25 it 0.05 for clusters 
at z ~ Ya, or JIq = 0.20 it 0.04. Our measurement of Q. is done completely within our data set 
which diminishes many of the possible selection effects to small corrections to the luminosity, 
filter bands, and shape of the luminosity function, since they are done identically for the field and 
cluster galaxies. The concern that relative evolution of field and cluster galaxies could lead to 



systematic errors is partially constrained here and is addressed elsewhere (Abraham et al. 1995 
and future papers). Cluster galaxies have largely ceased star formation, so if anything, are likely 
to be less luminous that field galaxies, meaning that the $7 given here is a mild overestimate. 

Complementary to the mean M^/L value is the average My/Ni, which is 3.5ib0.4x IQ^'^h^^ 
to M,.(0) = -18.5. and 5.7 ± 0.7 x 10^^ M^(0) = -19.5 (excluding the two binary clusters 
and MS1621-h26). The closure M/Nl values at these two luminosity levels give = 0.28 ± 0.08. 
Within the errors the M^/L and the M^/N^ (for two different limiting magnitudes) estimates 
for Q are consistent with a single value, but, because they are based on differing assumptions 
about the relative evolution of the cluster and field galaxy population, the sum squared errors 
in the two methods is used as an indicator of the systematic error resulting from the luminosity 
normalization, finding that it is 0.09. 
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Our self-contained technique eliminates many of the systematic errors of the M/L and M/N 

techniques for r2 estimation, but leaves open the issue whether cluster galaxies arc substantially 
more clustered than the cluster mass. This will be addressed with the same data set in a 
companion paper. 

We thank all participants of the CNOC cluster survey for assistance in obtaining and reducing 
these data. The Canadian Time Assignment Committee for the CFHT generously allocated 
substantial grants of observing time, and the CFHT organization provided the technical support 
which made these observations feasible. We gratefully acknowledge financial support from NSERC 
and NRC of Canada. 
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Fig. 1. — The redshift vs RA measured from the brightest cluster galaxy direction for the full fields 
of the 16 clusters. The redshift increases from the vertex of the diagram, and RA is converted to 
a transverse distance in physical co-ordinates. 

Fig. 2. — Redshift vs Declination, as in Figure 1. 

Fig. 3. — The redshift histograms in the vicinity of the clusters. The dashed lines indicate the 
selected cluster redshift range. 

Fig. 4. — The sky locations of the galaxies in the redshift range of the cluster. The circle area is 
proportional to the galaxy's apparent magnitude. Each box is 1000 arcseconds on a side. The same 
layout of the clusters as in Figures 1, 2 and 3q is used. North is to the top and East to the left. 

Fig. 5. — The corrected luminosities (top row) and normalized colors (bottom row) for cluster 
galaxies (left column) and field galaxies (right column). 

Fig. 6. — The redshift corrected constant overdensity luminosity plotted against the observed ai. 

Fig. 7. — The distribution of the deviations from the mean My/L (upper left, heavy lines for the 6 
well covered clusters), and the M^/L ratio as a function of total blue fraction, color, mean interior 
density, velocity dispersions and the redshift. 

Fig. 8. — The distribution of the deviations from the mean M^/Nl (upper left, heavy lines for 
the 6 well covered clusters) for Mj.{0) = —18.5. The M^/Nl ratio is plotted as a function of the 
quantiies the previous figure. 
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Table 1: Cluster Redshift and Field Limits 

Name z RA Dec zio Zhi 

h'^ Mpc 



A2390 


0. 


.2280 


6. 


.13 


1, 


.11 


0. 


.2175 


0.2395 


MS0016+16 


0, 


.5466 


2, 


.19 


1, 


.87 


0, 


.5300 


0.5587 


MS0302+16 


0, 


.4246 


1, 


.94 


1, 


.73 


0, 


.4160 


0.4300 


MS0440+02 


0, 


.1965 


3, 


.32 


1, 


.02 


0, 


.1910 


0.2013 


MS0451+02 


0, 


.2011 


4, 


.53 


1, 


.57 


0, 


.1919 


0.2120 


MS0451-03 


0, 


.5392 


2. 


.26 


1, 


.87 


0. 


.5240 


0.5530 


MS0839+29 


0, 


.1928 


3. 


.37 


0, 


.96 


0. 


.1820 


0.2024 


MS0906+11 


0, 


.1705 


1, 


.06 





.86 


0, 


.1567 


0.1838 


MS1006+12 


0, 


.2605 


1, 


.43 


1, 


.20 


0, 


.2530 


0.2680 


MS1008-12 


0, 


.3062 


1, 


.61 


1, 


.36 


0, 


.2934 


0.3144 


MS1224+20 


0, 


.3255 


1, 


.66 


1, 


.35 


0, 


.3180 


0.3310 


MS1231+15 


0, 


.2348 


1. 


.33 


3 


.16 


0. 


.2270 


0.2410 


MS1358+62 


0, 


.3290 


1, 


.72 


4, 


.27 


0. 


.3150 


0.3420 


MS1455+22 


0, 


.2568 


1, 


.47 


1, 


.21 


0, 


.2400 


0.2700 


MS1512+36 


0, 


.3726 


5, 


.51 


1, 


.70 


0, 


.3640 


0.3800 


MS1621+26 


0, 


.4274 


2, 


.04 


5, 


.15 


0, 


.4190 


0.4360 



Table 2: Velocity Dispersion Validation 
Name o-[/ai 1 s.d. interval 



A2390 


1. 


.07 


0. 


99 - 


- 1 


.11 


MS0016+16 


1, 


.07 


0. 


93 - 


- 1, 


.15 


MS0302+16 


0, 


.96 


0. 


,77 - 


- 1, 


.09 


MS0440+02 


0. 


.95 


0. 


60 ■ 


■ 1, 


.14 


MS0451+02 


0. 


.85 


0. 


,67- 


■ 0, 


.98 


MS0451-03 


0, 


.96 


0. 


,84 - 


- 1, 


.01 


MS0839+29 


0, 


.93 


0. 


,55 - 


- 1, 


.10 


MS0906+11 


0, 


.37 


0. 


22 - 


- 0, 


.49 


MS1006+12 


0, 


.93 


0. 


,74 - 


- 1, 


.00 


MS1008-12 


1. 


.03 


0. 


,90 ■ 


■ 1, 


.09 


MS1224+20 


0. 


.97 


0. 


,85 ■ 


■ 1, 


.04 


MS1231+15 


1, 


.04 


0. 


,95 - 


- 1, 


.10 


MS1358+62 


0, 


.96 


0. 


91 - 


- 1, 


.02 


MS1455+22 


0, 


.97 


0. 


,85 - 


- 1, 


.09 


MS1512+36 


1, 


.18 


1. 


00 - 


- 1, 


.32 


MS1621+26 


0. 


.92 


0. 


,83 ■ 


■ 1, 


.00 
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Table 3: CNOC Cluster Dynamical Parameters 



Name 


z 


N 






<7l 














Mpc 


kms 


-1 


Pc{z) 


A2390 


0.2280 


174 


1.972 


0.16 


1104 


63 


119 


MS0016+16 


0.5466 


47 


0.978 


0.12 


1234 


128 


360 


MS0302+16 


0.4246 


28 


0.502 


0.12 


639 


92 


442 


MS0440+02 


0.1965 


36 


1.224 


0.16 


606 


62 


99 


MS0451+02 


0.2011 


108 


1.371 


0.14 


1031 


73 


226 


MS0451-03 


0.5392 


51 


0.880 


0.10 


1371 


105 


555 


MS0839+29 


0.1928 


42 


0.558 


0.14 


756 


111 


745 


MS0906+11 


0.1705 


78 


0.584 


0.04 


1888 


117 


4426 


MS1006+12 


0.2605 


26 


0.632 


0.06 


906 


101 


738 


MS1008-12 


0.3062 


67 


0.575 


0.05 


1054 


107 


1113 


MS1224+20 


0.3255 


24 


0.478 


0.15 


802 


90 


903 


MS1231+15 


0.2348 


73 


0.944 


0.09 


640 


65 


173 


MS1358+62 


0.3290 


165 


1.472 


0.13 


934 


54 


128 


MS1455+22 


0.2568 


49 


0.621 


0.04 


1133 


151 


1203 


MS1512+36 


0.3726 


38 


1.019 


0.41 


690 


96 


136 


MS1621+26 


0.4274 


98 


1.428 


0.15 


793 


55 


84 



Table 4: CNOC Cluster and L 



Name 


p{rv) 




L 


M^/L 




M200 


L200 




Pc{z) 






hUQ 


/L© 






A2390 


119 


1.7x10^^ 


6.7x10^2 


248 


43 


1.3x10^5 


5.2x10^2 


MS0016+16 


360 


1.0x10^^ 


3.6x10^2 


288 


94 


1.3x10^^ 


4.8x10^2 


MS0302+16 


442 


1.4x10^'^ 


6.3x10^^ 


224 


108 


2.1x10^^ 


9.4x10" 


MS0440+02 


99 


3.1x10^"^ 


9.9x10" 


312 


87 


2.2x10^^ 


7.0x10" 


MS0451+02 


226 


1.0x10^^ 


2.8x10^2 


357 


71 


1.1x10^5 


3.0x10^2 


MS0451-03 


555 


1.1x10^^ 


2.9x10^2 


393 


102 


1.8x10^5 


4.8x10^2 


MS0839+29 


745 


2.2x10^4 


7.7x10" 


285 


83 


4.2x101-^ 


1.5x10^2 


MS0906+11 


4426 


1.4x10^5 


1.8x10^2 


800 


181 


6.6x10^5 


8.5x10^2 


MS1006+12 


738 


3.6x10^'^ 


1.2x10^2 


291 


92 


6.9x10^'' 


2.3x10^2 


MS1008-12 


1113 


4.4x10^^ 


2.0x10^2 


220 


60 


1.0x10^^ 


4.7x10^2 


MS1224+20 


903 


2.1x10^^ 


1.0x10^2 


212 


114 


4.5x10^^ 


2.1x10^2 


MS1231+15 


173 


2.7x10^^ 


1.5x10^2 


176 


43 


2.5x10^^ 


1.4x10^2 


MS1358+62 


128 


8.9x10^^ 


4.5x10^2 


197 


26 


7.1x10^^ 


3.6x10^2 


MS1455+22 


1203 


5.5x10^^ 


9.3x10" 


589 


201 


1.3x10^5 


2.3x10^2 


MS1512+36 


136 


3.3x10^^ 


1.4x10^2 


234 


156 


2.7x101"^ 


1.2x10^2 


MS1621+26 


84 


6.2x10^^ 


4.1x10^2 


151 


37 


4.0x10^^ 


2.7x10^2 
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Table 5: CNOC Cluster M/N 
Name M^/A^(-18.5) 



AzoyU 


z.oX iU 


Q 
O. 


MdUU1dH-1d 


O.oX iU 


z. 


MdUoUzH-1d 


Z.ZX iU 


Q 

o. 


i\ /Tcn/i A n 1 no 


O.oX iU 


y. 


MS0451+02 


3.8x10-'^^ 


5. 


MS0451-03 


7.6x10^2 


2. 


MS0839+29 


2.4x10^2 


7. 


MS0906+11 


9.6x10^2 


1. 


MS1006+12 


5.6x10^2 


1. 


MS1008-12 


2.6x10^2 


6. 


MS1224+20 


3.0x10^2 


1. 


MS1231+15 


1.8x10^2 


4. 


MS1358+62 


2.1x10^2 


2. 


MS1455+22 


5.7x10^2 


1. 


MS1512+36 


2.2x10^2 


1. 


MS1621+26 


2.0x10^2 


4. 





M/Ar(-19.5) 


Cm/at 


1,-1 


Mq 




6x10^^ 


5.7x10^^ 


1.1x10^^ 


0x10 


o > . "1 r\'\ 2 

5.8x10 


2.0x10 


6x10^^ 


3.4x10 


2.1x10 


2x10^^ 


6.8x10^^ 


2.3x10^2 


9x 10^^ 


7.4x10^^ 


1.7x10 


n 1 nl 2 

0x10^^ 


7.6x10^^ 


o 1 1 nl 2 

2.1x10 


6x10" 


7.7x10^2 


3.2x10^2 


9x10^2 


2.2x10^3 


6.1x10^2 


8x10^2 


5.6x10^2 


1.8x10^2 


9x10" 


5.0x10^2 


1.7x10^2 


6x10^2 


4.6x10^2 


2.1x10^2 


5x10" 


4.1x10^2 


1.3x10^2 


6x10" 


4.2x10^2 


6.2x10" 


7x10^2 


1.4x10^3 


5.8x10^2 


5x10^2 


6.7x10^2 


4.4x10^2 


1x10" 


2.6x10^2 


4.7x10" 



Table 6: Average My/L and M^/Nl varying the Absolute Magnitude Limit 



M^(0) 


M,JL 










hUQ 


/L© 




Mo 


-18.0 


310 


113 


3.3 X 10^2 


1.7 X 10^2 


-18.5 


283 


98 


3.5 X 10^2 


1.6 X 10^2 


-19.0 


247 


80 


3.9 X 10^2 


1.7 X 10^2 


-19.5 


241 


90 


5.7 X 10^2 


2.4 X 10^2 



Table 7: The Closure M/L and M/Nl Values varying the Absolute Magnitude Limit 

^r(O) Pc/j{> L) a M/L Pcl4>{> L) Um/N 

/iMq/L© h-^M.Q 

-18.0 1265 149 1.22 x lO^^ 2.I x 10^2 

-18.5 1163 130 1.32 x 10^^ 1.9 x 10^2 

-19.0 1073 104 1.60 x 10^^ 1.5 x 10^2 

-19.5 1007 103 2.37 x 10^^ 2.6 x 10^2 



